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A modification of the spiked harmonic oscillator is studied in the case for which the per- 
turbation potential contains both an inverse power and a linear term. It is then possible 
to evaluate trial functions by solving an integral equation due to the occurrence of the 
linear term. The general form of such integral equation is obtained by using a Green- 
function method, and adding a modified Bessel function of second kind which solves an 
homogeneous problem with Dirichlet boundary condition at the origin. 
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1. INTRODUCTION 



In the literature on quantum-mechanical problems, an important role is played by the 
"spiked" harmonic oscillator in three spatial dimensions. This is a system where the radial 
part ip{r) of the wave function is ruled by the Hamiltonian operator [1-4] 

7 TI ,n d2 2 1(1 + 1) A 

in the sense that H(a, A) acts on cp(r) = rvjj(r). In the s-wave case (i.e. when the angular 
momentum quantum number I vanishes) H(a, A) reduces to 

If (a, A) = -^ + r 2 + A = F (a, 0) + AV, (1.2) 

where H(a, 0) is formally equal to the simple harmonic oscillator Hamiltonian, r belongs 
to the interval [0, oo] and V = r~ a . For any fixed value of A, the potential term in (1.2) 
diverges as r — > in such a way that the operator if (a, A) acts on wave functions which 
vanish at the origin: 

p(0) = 0. (1.3) 

More precisely, the imposition of Eq. (1.3) is necessary since not all functions in the 
domain of H(a,0) are in the domain of V [2]. Thus, when A — > and a is fixed, the 
operator H(a, A) converges to an operator formally equal to the unperturbed operator 
Hq = —-^2 + 1 but supplemented by the Dirichlet boundary condition (1.3) for all 
functions in its domain. This means that the unperturbed operator Hq, for which the 
boundary condition (1.3) is not necessary, differs from the limiting operator H(a,0), for 
which Eq. (1.3) is instead necessary to characterize the domain. 

The full potential in (1.1) or (1.2) inherits the name "spiked" from a pronounced peak 
near the origin for A > 0, and its consideration is suggested by many concrete problems in 
chemical, nuclear and particle physics. Some important results within this framework are 
as follows. 

(i) Development of singular perturbation theory, with application to the small- A expansion 
of the ground-state energy [1] . 

(ii) A variational method has been successfully applied to a large-coupling perturbative 
calculation of the ground-state energy [2] . 
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(iii) Weak-coupling perturbative analysis [3], and its relation with the strong coupling 
regime for a in the neighbourhood of a = 2. 



(iv) A non-perturbative but absolutely convergent algorithm for the evaluation of eigen- 
functions [4]. 

Note now that for a > 2 the Hamiltonian operators (1.1) or (1.2) lead to a non- 
Fuchsian singularity [5] at r = of the stationary Schrodinger equation, because their 
potential term has a pole of order > 2 therein. Section 2 outlines the method developed 
by Harrell [1] for dealing with such singularities in a perturbative analysis. A non-trivial 
extension is studied in Sec. 3, i.e. a model where the perturbation potential contains both 
an inverse power and a linear term. Concluding remarks are presented in Sec. 4. 

2. THE HARRELL METHOD 

The method developed by Harrell relies on the choice of suitable trial functions for self- 
adjoint operators (i.e. normalized vectors in their domain) and on the following lemma 



Lemma. If ip\ is a trial function for the self-adjoint operator T + XT' , where both T and 
T' are self-adjoint and E(0) is an isolated, non-degenerate stable eigenvalue of T, and E(X) 

is a continuous function such that the scalar product [T + XT' — E(\)]ip\j tends to 

as A tends to 0, and 



In a one-dimensional example, hereafter chosen for simplicity, let the full Hamiltonian 
be H + XV, where 




|| [ T + AT' - E(XM X \\ = o (V(Va, [T + XT' - E(X)]^ x )"j 



(2.1) 



then the eigenvalue of T + XT' which converges to E(0) is 



E(X) =(Va, [T + AT']Va) + O (||[T + AT' - E(X)]ip\\\ 2 ^ 



(2.2) 



#0 = 



dx 2 




(2.3) 



and 



V(x) 



= x 



— a 



(2.4) 
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If a = 4, on denoting by Ui the unperturbed eigenstates, a trial function can be chosen in 
the form [1] 

i{;(x) = W(x; X)ui(x), (2.5) 



where (inspired by the JWKB approximation) 



W(x; A) = iV(A)exp (- y/W{£)d^ = iV(A)exp 



x 



(2.6) 



with N(X) a normalization factor that approaches 1 as A tends to 0. One then finds the 
formula [1] 



[H -Ei + \V]W Ui = 2W^- (±- - 



Ui. 



(2.7) 



By virtue of Eq. (2.7) and of the Lemma at the beginning of this section one obtains 
the following formula for the energy eigenvalues of the spiked harmonic oscillator in one 
dimension with a = 4: 

Ei(X) = Ei(0) + 2 (ui,x~ 2 Q - ^ u^j V\ + 0(A). (2.8) 



When a^4, one can use the identity (hereafter W a replaces W) 

d 2 



dx 2 



(w aUi ) = w a u'l + 2W>J + wZm, 



(2.9) 



which implies [1] 



[Hq - Ei + XV]W aUi = 



d 2 W n dW n d 



dx 2 



dx dx 



+ xvw a 



Ui. 



(2.10) 



The idea is now to choose W a in such a way that the action of [Hq — E { + XV] on W a Ui 
involves again the action of the operator (^ — ^) on Ui (see Eq. (2.7)). For this purpose, 
Harrell imposed the differential equation 



d 2 2d 
+ 



X 



dx 2 x dx x a 



W a = 0, 



(2.11) 



so that Eq. (2.10) reduces indeed to (cf. Eq. (2.7)) 



[H -E i + XV]W a u i = 2- 



dW n (\ d 



dx V x dx 



Ui 



(2.12) 



On denning 



the solution of Eq. (2.11) can be expressed in the form 

W a {x-X) = ^^-x-^K^vVlx-^). (2.14) 

The energy eigenvalues of the spiked oscillator in one dimension with a > 4 are then found 
to be [1] 

, r (! ~ v ),2v (.. ^-if 1 d \ „. A x» , r^^2u, 



/•-/I A) = /•:,(()) + 2 r | [ ^V " [u^x 1 \- .. _)„,) A" • O(A^). (2. 1 5) 



3. EXTENSION TO OTHER SINGULAR POTENTIALS 

When the perturbation potential is not an inverse power, Eq. (2.11) is replaced by an 
equation which cannot generally be solved. Harrell has however shown that, if V is bounded 
away from x = and lies in between x~ a and x~@ : with < a < (3, then the effect of V 
on the eigenvalues is not essentially different [1] . 

It therefore appears of interest to study cases in which V is not (a pure) inverse power, 
but does not obey the restrictions considered in Sec. 5 of Ref. [1] . For this purpose, we 
assume that (cf. (2.4)) 

V(x) = x~ a + kx. (3.1) 

By doing so we study a model that reduces to the spiked oscillator in the neighbourhood 
of the origin, whereas at large x it approaches an oscillator perturbed by a "Starkdike" 
term. The two terms in (3.1) are separately well studied in the literature, so that their 
joint effect provides a well motivated departure from the scheme studied in Sec. 5 of Ref. 

[!]• 

The exact solution of the counterpart of Eq. (2.11) when (3.1) holds is not available 
in the literature to our knowledge. For this purpose we study an integral equation, whose 
construction is as follows. On assuming the validity of Eq. (3.1) for the perturbation 
potential, Eq. (2.11) is replaced by the inhomogeneous equation 

LW a (x;X) = f a (x;X), (3.2) 
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where 

d 2 2 d A 
L= — + --: , 3.3 

ax z x ax x a 

fa(x. A) = \KxW a (x; A). (3.4) 

Since V(x) approaches x~ a as x — > 0, we require again the boundary condition studied in 
Ref. [1], i.e. 

W a (0) = 0. (3.5) 

We are here following an approach similar to the one leading to an integral equation 
for scattering problems, where the right-hand side of the stationary Schrodinger equation 
involves again the unknown function: 

ZT77 

(/\ + k 2 )Mx) = - T v(x)MS). 

n 

In our problem, following Ref. [1], we study one-dimensional Schrodinger operators with 
Dirichlet boundary conditions at only on the space L 2 ([0, oo)), to remove the degeneracy 
resulting from the decoupling of the two half-lines (— oo, 0) and (0, oo) (only odd eigen- 
functions of the ordinary harmonic oscillator obey the Dirichlet condition at 0). Moreover, 
to be able to use all known results on one-dimensional boundary- value problems on closed 
intervals of the real line, we first work on the interval (0, b) and then take the limit as 
b — > oo. More precisely, we start from a problem whose Green function G(x,£; A) satisfies 
the equation 

LG = for xe (0, oo) and f e (0, oo), (3.6) 

the boundary condition 

G(0,fcA) = 0, (3.7) 

the summability condition (in that the integral of G over a closed interval of values of A 
yields a square-integrable function of x) 

G(x,£;A)e4 c) (0,oo), (3.8) 



and the continuity conditions 



lim G(x,£;A)= lim G(x,£;A), (3.9) 



lim ™ _ li m ^ = 1. (3.10) 

OX as->-£- OX 
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As is shown on page 442 of Ref. [5], G(x,£; A) is recovered on studying first the Green 
function G b (x, £; A) for a regular problem on the interval (0, 6), and then taking the limit 
as b — > oo, i.e. 

hm G b (x,e;A) = G(x,e;A). (3.11) 

Such a relation holds independently of the boundary condition imposed on G b (x,£; A) at 
x = b [5]. With this understanding, the full solution of the inhomogeneous equation (3.2) 
with boundary condition (3.5) is given by (7 being a constant) 

W a (x;X) = 7 2 p^ x~*K u (2vV\x~&) 

+ hm Xk [ G b (x,Z;\)ZW a (b\)dt, (3.12) 

where the first term on the right-hand side of (3.12) is the regular solution (2.14) of 
the homogeneous equation LW a = 0. The Green function G b (x,£;\) has to obey the 
differential equation [5] 

LG b = for x E (0,f) and x E (£,&), (3.13) 
the homogeneous boundary conditions [5] 

G b ((U;A) = 0, (3.14) 

G b (b,£;\) = 0, (3.15) 

and the continuity conditions [5] 

lim G b (x,£;\)= lim G b (x,£;\), (3.16) 

lim — lim = 1. (3-17) 

x^£+ OX x^>£- OX 

To obtain the explicit form of G b (x, £; A) one has to consider a non-trivial solution uq(x; A) 
of the homogeneous equation Lu = satisfying u(0) = 0, and a non-trivial solution u b (x; A) 
of Lu = satisfying u(b) = 0. By virtue of (3.13)-(3.15) one then finds [5] 

G b (x, £; A) = X)u (x; A) if xG (0, £), (3.18) 

G 6 (x,£;A) = B(£;A)u 6 (x;A) if iG^ft), (3.19) 
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where uq and u b are independent. The continuity conditions (3.16) and (3.17) imply that 
A(£; A) and £?(£; A) are obtained by solving the inhomogeneous system 



A(£;A)«o(£;A)-£?(£;AH(£;A) = 0, 
B&X)u b &X)-A&X)u' &X) = l, 



which yields 



A) = 



Ub(& X) 



(3.20) 
(3.21) 

(3.22) 
(3.23) 



where O is the Wronskian of uq and u b . We now recall the Abel formula for fi, according 
to which [5] 

Q(u , u h - x; A) = C(X)e- v{x \ (3.24) 



where C(A) is a constant and, for the operator L in (3.3), v is a particular solution of the 
equation [5] 

dv _ 2 
dx x ' 



i.e. 



which implies that 



v(x) = log a; 2 , 
C 



(3.25) 
(3.26) 



(3.27) 



Thus, on defining as usual x< = min(x,^),x> = max(x,^), Eqs. (3.18), (3.19), (3.22), 
(3.23) and (3.27) imply that the Green function is expressed by 



C(X) 



w (^<; X)u b (x > ; A). 



(3.28) 



The integral equation (3.12) for W a becomes therefore 



W a {x;X) = 1 —-—x--*K v (2v^f\x-^) 



,. Ak 
hm — 



r(iz) 

u b (x;X) [ «o(^;A)e 3 W Q (e;A)^ 
Jo 



+u (x;X) / u b (C;X)eW ce (C;X)di 



(3.29) 
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In Eq. (3.29), uq(x; A) and Ub(x; A) can be chosen to be of the form 



u (x; A) = Cq{u)x ^ K l/ (2uy/\x 2 -) 



(3.30) 



Ub(x; A) = Cb(v)x 2/ I/ (2z/v / Ax 



(3.31) 



in agreement with the homogeneous Dirichlet conditions at and at b, respectively. Before 
taking the limit as b — > oo in Eq. (3.29) we can now regard the perturbation parameter A 
as an eigenvalue. We are therefore studying, for finite b, a Fredholm integral equation of 
second kind, whose general form is (the parameter a vanishes in our problem) 



If A is an eigenvalue, a necessary and sufficient condition for the existence of solutions of 
Eq. (3.32) is that, for any solution \ °f the equation 



In our problem, / is the first term on the right-hand side of Eq. (3.29), the kernel K is 
given by Gb in (3.28), and Eq. (3.34) provides a powerful operational criterion. Indeed, 
one might try to use directly the theory of integral equations [6] on the interval (0, oo) 
instead of the limiting procedure in Eq. (3.29), but the necessary standard of rigour goes 
beyond our present capabilities. 

4. CONCLUDING REMARKS 

In the present letter we have exploited the fact that if a spiked harmonic oscillator is 
modified by the addition of a linear term, the full perturbation potential may be seen as 
consisting of an inverse power plus a term linear in the independent variable. It is then 
possible to evaluate the function W(x; A) occurring in the trial function (2.5) by solving 
the inhomogeneous equation (3.2), which leads to the integral equation (3.29). This is 
involved, but leads in principle to a complete calculational scheme. Note that, if one 
tries to combine the term x 2 in the unperturbed Hamiltonian with the linear term in the 




(3.32) 




(3.33) 



the known term f(s) should satisfy the condition [6] 




(3.34) 
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perturbation potential (3.1), one eventually moves the singular point away from the origin, 
whereas the spiked oscillator is (normally) studied by looking at the singular point at the 
origin. 

It would be rather interesting, as a subject for further research, to consider suitable 
changes of independent variable in the investigation of non-Fuchsian singularities. For 
example, given the Hamiltonian operator 



^ _ d b a 



(4.1) 



if one defines the new independent variable (cf. page 971 of Ref. [7]) 



p = r 



(4.2) 



for a suitable parameter 7, the stationary Schrodinger equation becomes 



(f ( 1\ 1 d 1/2 (2-p) X 

^+[ 1 -t )z- + — ( e P-> -<>P ~< 



7 / p dp 7 2 p 2 



<p(p) = 0. 



(4.3) 



By construction, the larger is 7, the more Eq. (4.3) tends to its Fuchsian limit 



d 2 f l 1^ 1 d | (E-a-b) 1 



dp 



+ 



l) pdp 



r 



<p(p) = 0, 



(4.4) 



for all values of p. This remark can be made precise by defining 



1 

£ = -, 

7 



F(e) =e 2 {Ep 2s -ap^ £ -b), 
and considering the asymptotic expansion at small e (and hence large 7) 



(4.5) 
(4.6) 



F(e) ~e 2 E - a - b + e(2E - (2 - p)a) log p + 0(e 2 ) 



(4.7) 



The first "non-Fuchsian correction" of the limiting equation (4.4) is therefore 



^ + (1 " e)~ + j 2 [e 2 ((E - a - b) + e(2E - (2 - p)a) logp) } <p(p) = 0. (4.8) 
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Interestingly, logarithmic terms in the potential can be therefore seen to result from a 
sequence of approximations relating Eq. (4.3) to its Fuchsian limit (4.4). Moreover, all 
equations with Fuchsian singularities like Eq. (4.4) might be seen as non-trivial limits of 
stationary Schrodinger equations with non-Fuchsian singular points. It remains to be seen 
whether such properties can be useful in the investigation of the topics discussed in the 
previous sections. 

Another topic for further research is the Schrodinger equation for perturbed stationary 
states of an isotropic oscillator in three dimensions written in the form 



— + k 2 

dr 2 



2 2 

/i r 



1(1 + 1) 



Sir) 



<p(r) = 0, 



(4.9) 



where, having set (here fx = ^ 



V(r) = -^U(r)=/j 2 r 2 + S(r), 

1% 



(4.10) 



the function S represents the "singular" part of the potential according to our terminology. 
We look for exact solutions of Eq. (4.9) which can be written as 



tp(r) = A(r)e s(r) e-*S-. 



(4.11) 



The second exponential in (4.11) takes into account that, at large r, the term A 2 r 2 domi- 
nates over all other terms in the potential (including, of course, ^ s ), and has not been 
absorbed into B(r) for later convenience. It is worth stressing that Eq. (4.11) is not a 
JWKB ansatz but rather a convenient factorization of the exact solution of Eq. (4.9). 
We determine B(r) from a non-linear equation by straightforward integration (see below), 
while the corresponding second-order equation for A is rather involved. 
Indeed, insertion of (4.11) into Eq. (4.9) leads to 



+ 2(B' -vr)^+\k 2 -v- - 2»rB' + B" + B' 2 - S(r) 



A(r) = 0. 



(4.12) 

To avoid having coefficients of this equation which depend in a non-linear way on B we 
choose the function B so that 



B' 2 - S(r) = 0, 



(4.13) 
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which implies (up to a sign, here implicitly absorbed into the square root) 



B(r) = J ^fS{?)dr. (4.14) 
Hence one finds the following second-order equation for the function A: 

{ £ + 2( ^ - H ± + [k 2 - „ - ^±12 - 2/uVS + ^=J J A(r) = 0. (4.15) 

It should be stressed that the step leading to Eq. (4.13) is legitimate but not compelling. 
For each choice of B(r) there will be a different equation for A(r), but in such a way that 
<p(r) remains the same (see (4.11)). Unfortunately, Eq. (4.15) remains too difficult, as far 
as we can see. 
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